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A significant improvement to the development of CFD-based unsteady aerodynamic 
reduced-order models (ROMs) is presented. This improvement involves the simultaneous 
excitation of the structural modes of the CFD-based unsteady aerodynamic system that 
enables the computation of the unsteady aerodynamic state-space model using a single 
CFD execution, independent of the number of structural modes. Four different types of 
inputs are presented that can be used for the simultaneous excitation of the structural 
modes. Results are presented for a flexible, supersonic semi-span configuration using the 
CFL3Dv6.4 code. 


Introduction 

E ARLY mathematical models of unsteady aerody- 
namic response capitalized on the efficiency and 
power of superposition of scaled and time-shifted fun- 
damental responses, also known as convolution. Clas- 
sical models of two-dimensional airfoils in incompress- 
ible flow 1 include Wagner’s function 2 (response to a 
unit step variation in angle of attack), Kussner’s func- 
tion 3 (response to a sharp-edged gust), Theodorsen’s 
function 4 (frequency response to sinusoidal pitching 
motion), and Sear’s function 5 (frequency response to 
a sinusoidal gust). As geometric complexity increased 
from airfoils to wings to complete configurations, the 
analytical derivation of these types of response func- 
tions became impractical and the numerical computa- 
tion of linear unsteady aerodynamic responses, in the 
frequency domain, became the method of choice . 6 

When geometry- and flow-dependent nonlinear aero- 
dynamic effects became significant, appropriate non- 
linear aerodynamic equations were solved using time- 
integration techniques. Coupling the nonlinear aero- 
dynamic equations with a linear structural model pro- 
vides a direct simulation of aeroelastic phenomena. 
This direct simulation approach for solving nonlinear 
aeroelastic problems has yielded a very powerful sim- 
ulation capability with two primary challenges. The 
first challenge is the associated computational cost 
of this simulation, which increases with the fidelity 
of the nonlinear aerodynamic equations to be solved. 
Computational cost may be reduced via the imple- 
mentation of parallel processing techniques, advanced 
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algorithms, and improved computer hardware process- 
ing speeds. The second, more serious, challenge is that 
the information generated by these simulations cannot 
be used effectively within a preliminary design envi- 
ronment. Any attempt to incorporate the output of 
these aeroelastic simulations within a design environ- 
ment inevitably becomes design by trial-and-error. As 
a result, the integration of traditional, computational 
aeroelastic simulations into preliminary design activi- 
ties involving disciplines such as aeroelasticity, aeroser- 
voelasticity (ASE), and optimization is, at present, a 
costly and impractical venture. 

The goal behind the development of reduced-order 
models (ROMs) is aimed precisely at addressing these 
two challenges. Development of a ROM entails the 
development of a simplified mathematical model that 
captures the dominant dynamics of the original sys- 
tem. This alternative mathematical representation of 
the original system is, by design, in a mathematical 
form suitable for use in a multidisciplinary, prelimi- 
nary design environment. As a result, interconnection 
of the ROM with other disciplines is possible, thereby 
addressing the second challenge. The simplicity of 
the ROM yields significant improvements in compu- 
tational efficiency as compared to the original system, 
thereby addressing the first challenge. 

At present, the development of CFD-based ROMs 
is an area of active research at several government, 
industry, and academic institutions . 7-11 Development 
of ROMs based on the Volterra theory is one of sev- 
eral ROM methods currently under development . 12-16 
Reduced-order models based on the Volterra theory 
have been applied successfully to Euler and Navier- 
Stokes models of nonlinear unsteady aerodynamic and 
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aeroelastic systems. Volterra-based ROMs are based 
on the creation of linearized and nonlinear unsteady 
aerodynamic impulse responses that are then used in a 
convolution scheme to provide the linearized and non- 
linear responses of the system to arbitrary inputs. In 
this setting, the linearized and nonlinear impulse re- 
sponses may be considered to be the ROMs of the 
particular nonlinear system under investigation. Al- 
ternatively, upon transformation of the linearized and 
nonlinear impulse responses into state-space form, the 
state-space models generated may also be considered 
ROMs. 

Various inputs can be used in the time domain 
(CFD code) to generate generalized aerodynamic 
forces (GAFs) in the frequency domain in order to 
perform standard, frequency-domain aeroelastic anal- 
yses. But if time-domain aeroservoelastic (ASE) anal- 
yses are desired, the frequency-domain GAFs are 
transformed back into the time domain using tradi- 
tional rational function approximation (RFA) tech- 
niques. These techniques include, for example, the 
well-known Rogers approximation 17 and the Minimum 
State technique. 18 The RFA techniques transform 
frequency-domain GAFs into state-space (time do- 
main) models amenable for use with modern control 
theory and optimization. The process just described 
transforms time-domain information (CFD results) 
into frequency-domain information only to have the 
frequency-domain information transformed back into 
the time domain. 

Gupta et al 19 and Cowan et al 20,21 applied a set 
of flight testing inputs to an unsteady CFD code and 
used the information to create a linear autoregressive 
moving average (ARMA) model that was transformed 
into state-space form. Although this technique is ap- 
plied entirely within the time domain, the shape of the 
inputs applied to the CFD code requires tailoring in 
order to excite a specific frequency range, resulting in 
an iterative process. In a similar vein, Rodrigues 22 de- 
veloped a state-space model for an airfoil in transonic 
flow using a transonic small-disturbance algorithm. 

In the present paper, a direct approach for effi- 
ciently generating linearized unsteady aerodynamic 
state-space models is presented. Although the present 
application of the method deals with linearized re- 
sponses based on linearized impulse responses (lin- 
earized Volterra kernels), the method can be formally 
extended to address nonlinear aeroelastic phenomena 
via the use of nonlinear impulse responses (nonlinear 
Volterra kernels). 

Silva and Bartels 8 introduced the development of 
linearized, unsteady aerodynamic state-space models 
for prediction of flutter and aeroelastic response using 
the parallelized, aeroelastic capability of the CFL3Dv6 
code. The results presented provided an important 


validation of the various phases of the ROM devel- 
opment process. The Eigensystem Realization Algo- 
rithm (ERA), 23 which transforms an impulse response 
(one form of ROM) into state-space form (another 
form of ROM), was applied for the development of 
the aerodynamic state-space models. The ERA is 
part of the SOCIT (System/Observer/Controller Iden- 
tification Toolbox). Flutter results for the AGARD 
445.6 Aeroelastic Wing using the CFL3Dv6 code were 
presented as well, including computational costs. Un- 
steady aerodynamic state-space models were gener- 
ated and coupled with a structural model within a 
MATLAB/SIMULINK 24 environment for rapid calcu- 
lation of aeroelastic responses including the prediction 
of flutter. Aeroelastic responses computed directly 
using the CFL3Dv6 code showed excellent compari- 
son with the aeroelastic responses computed using the 
CFD-basecl ROM. 

Previously, the aerodynamic impulse responses that 
were used to generate the aerodynamic state-space 
models were computed using CFL3Dv6.0 via the exci- 
tation of one mode at a time. For a four-mode system, 
these computations are not very expensive. However, 
for more realistic cases where the number of modes 
can be an order of magnitude or more larger, the 
one-mode-at-a-time method becomes impractical. To- 
wards the solution of this problem, Raveh, 11 Kim 25 
and Kim et al 10 developed methods that enable the si- 
multaneous excitation of the structural modes, greatly 
reducing the cost of identifying the aerodynamic im- 
pulse responses from the CFD code. Raveh’s method 
consists of using filtered white Gaussian noise as si- 
multaneous excitation to the structural modes. Raveh 
then presents three methods for generating ROMs: a 
frequency-domain method, an ARMA method, and a 
state-space method. Kim’s method consists of using 
staggered step inputs, one per mode, and then ap- 
plying several techniques for recovering the individual 
responses from this type of simultaneous excitation 
leading to the creation of state-space models. 

This paper introduces four different types of ex- 
citations that can be used to simultaneously excite 
the structural modes while enabling the recovery of 
the individual responses. This will enable the com- 
putation of aerodynamic impulse responses for any 
number of structural modes using a single CFD ex- 
ecution. Reduced-order models generated using the 
one-mode-at-a-time method are compared to ROMs 
using the simultaneous excitation inputs. Recent ad- 
ditional enhancements to the development of unsteady 
aerodynamic and aeroelastic ROMs, which includes 
the use of ROMs for static aeroelastic responses at 
matched-point atmospheric conditions, are presented 
in a separate paper. 26 
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Description of CFD and System 
Identification Methods 

The following subsections describe the parallelized, 
aeroelastic version of the CFL3Dv6.4 code, the phases 
of the original and improved ROM development pro- 
cesses, and a description of the four types of functions 
used for simultaneous excitation of the CFD unsteady 
aerodynamic system. 

CFL3Dv6.4 Code 

The computer code used in this study is the 
CFL3Dv6 code, which solves the three-dimensional, 
thin-layer, Reynolds averaged Navier-Stokes equations 
with an upwind finite volume formulation. 27-29 The 
code uses third-order upwind-biased spatial differenc- 
ing for the inviscid terms with flux limiting in the 
presence of shocks. Either flux-difference splitting or 
flux- vector splitting is available. The flux-difference 
splitting method of Roe 30 is employed in the present 
computations to obtain fluxes at cell faces. There are 
two types of time discretization available in the code. 
The first-order backward time differencing is used for 
steady calculations while the second-order backward 
time differencing with subiterations is used for static 
and dynamic aeroelastic calculations. Furthermore, 
grid sequencing for steady state and multigrid and lo- 
cal pseudo-time stepping for time marching solutions 
are employed. 

One of the important features of the CFL3D code is 
its capability of solving multiple zone grids with one- 
to-one connectivity. Spatial accuracy is maintained 
at zone boundaries, although subiterative updating of 
boundary information is required. Coarse-grained par- 
allelization using the Message Passing Interface (MPI) 
protocol can be utilized in multiblock computations by 
solving one or more blocks per processor. When there 
are more blocks than processors, optimal performance 
is achieved by allocating an equal number of blocks to 
each processor. As a result, the time required for a 
CFD-basecl aeroelastic computation can be dramati- 
cally reduced. 

In this paper, multiblock MPI parallel aeroelastic 
computations, including flutter, for a flexible, semi- 
span supersonic configuration are performed using 28 
flowfield blocks. In order to achieve an optimal di- 
vision of grid points, it is necessary to place flow 
field block boundaries near a moving solid surface (the 
wing). The multiblock boundary and interior move- 
ment scheme allows the user to place block boundaries 
near surfaces as necessary for optimal parallelization. 
Boundaries interior to the fluid domain near a sur- 
face respond to the local surface motion. As the wing 
moves, block boundaries move to maintain integrity of 
block interfaces and the airfoil surface. 

Because the CFD and computational structural me- 


chanics (CSM) meshes usually do not match at the 
interface, CFD/CSM coupling requires a surface spline 
interpolation between the two domains. The interpo- 
lation of CSM mode shapes to CFD surface grid points 
is done as a preprocessing step. Modal deflections at 
all CFD surface grids are first generated. Modal data 
at these points are then segmented based on the split- 
ting of the flow field blocks. Mode shape displacements 
located at CFD surface grid points of each segment are 
used in the integration of the generalized modal forces 
and in the computation of the deflection of the de- 
formed surface. The final surface deformation at each 
time step is a linear superposition of all the modal 
deflections. 

System Identification Method 

In structural dynamics, the realization of discrete- 
time state-space models that describe the modal dy- 
namics of a structure has been enabled by the de- 
velopment of algorithms such as the Eigensystem 
Realization Algorithm (ERA) 23 and the Observer 
Kalman Identification (OKID) 31 Algorithm. These 
algorithms perform state-space realizations by us- 
ing the Markov parameters (discrete-time impulse 
responses) of the systems of interest. These algo- 
rithms have been combined into one package known as 
the System/Observer/Controller Identification Tool- 
box (SOCIT) 32 developed at NASA Langley Research 
Center. 

There are several algorithms within the SOCIT that 
are used for the development of unsteady aerody- 
namic discrete-time state-space models. The PULSE 
algorithm is used to extract individual input/output 
impulse responses from simultaneous input/output re- 
sponses. For a four-input/four-output system, simul- 
taneous excitation of all four inputs 33 yields four out- 
put responses. The PULSE algorithm is used to ex- 
tract the individual sixteen (four times four) impulse 
responses that associate the response in one of the out- 
puts due to one of the inputs. Details of the PULSE 
algorithm are provided in the references. Once the in- 
dividual sixteen impulse responses are available, they 
are then processed via the Eigensystem Realization 
Algorithm (ERA) in order to transform the sixteen 
individual impulse responses into a four-input/four- 
output, discrete-time, state-space model. A brief sum- 
mary of the basis of this algorithm follows. 

A finite dimensional, discrete-time, linear, time- 
invariant dynamical system has the state-variable 
equations 

x(k + 1) = Ax{k) + Bu{k ) (1) 

y(k) = Cx(k) + Du(k ) (2) 

where x is an n-dimensional state vector, u an m- 
dimensional control input, and y a p-dimensional out- 
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put or measurement vector with k being the discrete 
time index. The transition matrix, A, characterizes 
the dynamics of the system. The goal of system real- 
ization is to generate constant matrices (A, B, C) such 
that the output responses of a given system due to a 
particular set of inputs is reproduced by the discrete- 
time state-space system described above. 

For the system of Eqs. (1) and (2), the time-domain 
values of the systems discrete-time impulse response 
are also known as the Markov parameters and are de- 
fined as 

Y{k) = CA k ~ 1 B (3) 

with B an (n x m) matrix and C a (p x n) matrix. 
The ERA algorithm begins by defining the generalized 
Hankel matrix consisting of the discrete-time impulse 
responses for all input/output combinations. The al- 
gorithm then uses the singular value decomposition 
(SVD) to compute the A, B, and C matrices. 

In this fashion, the ERA is applied to unsteady 
aerodynamic impulse responses to construct unsteady 
aerodynamic state-space models. 

Original ROM Development Processes 

A CFD-based aeroelastic system can be viewed as 
the coupling of a nonlinear unsteady aerodynamic sys- 
tem (flow solver) with a structural system as depicted 
in Figure 1. The present study focuses on the devel- 
opment of a linearized unsteady aerodynamic ROM 
(in state-space form), using the general procedure de- 
picted in Figure 2, that is then coupled to a structural 
model (also in state-space form) for aeroelastic analy- 
ses. For the discussions that follow, the term ROM 
will refer to the unsteady aerodynamic state-space 
model. When the unsteady aerodynamic state-space 
model (ROM) is connected to a state-space model of 
the structure within the SIMULINK environment, this 
system is often also referred to as a ROM. However, 
to avoid confusion, the SIMULINK aeroelastic system 
will be referred to as the aeroelastic simulation ROM. 

An outline of the original (one-mode-at-a-time) 
ROM development process 8 is presented as back- 
ground for the new enhancements. The original ROM 
development process is as follows: 

1. Implementation of impulse/step response tech- 
nique into aeroelastic CFD code; 

2. Computation of impulse/step responses for each 
mode, one mode at a time, of an aeroelastic system 
using the aeroelastic CFD code; these responses are 
computed about a static aeroelastic solution (or a 
given dynamic pressure); 

3. Impulse responses generated in Step 2 are trans- 
formed into an unsteady aerodynamic state-space sys- 
tem using the ERA (within SOCIT); 

4. Evaluation/ validation of the state-space models 
generated in Step 3 via comparison with CFD results 



Fig. 1 Coupling of structure and aerodynamics 
within an aeroelastic CFD code. 


l 

' LINEAR 
' STRUCTURE 



Fig. 2 Identification of generalized aerodynamic 
forces (GAFs). 

(i.e. , ROM results vs. full CFD solution results); 

In the original ROM process, since each mode is be- 
ing excited individually, the response in each output 
due to a particular input (e.g., the impulse response 
in output 2 due to input 1) is generated almost di- 
rectly. If the function being used to excite the system 
is an impulse (or a unit pulse for discrete-time sys- 
tems), then the output from the CFD solution will 
consist of the impulse response of each output due 
to that one input. If the function being used to ex- 
cite the system is not an impulse function (such as a 
step input or a random input), then the impulse re- 
sponse needs to be extracted from the input/output 
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data. One method for extraction is deconvolution 8 
but there are other methods that can be used. For the 
original ROM process, since each input to a system 
is being excited individually, the exact nature of the 
input function is defined, to a certain extent, based 
on user preference (frequency range of interest, ease 
of implementation into a CFD solver, etc.). The step 
input has emerged as a convenient input for the orig- 
inal ROM process due to its ease of implementation 
into a CFD solver and to the wide range of frequency 
excitation that it can generate. 

The primary issue, then, with the original ROM pro- 
cess was the identification of the unsteady aerody- 
namic impulse responses one mode at a time (Step 2). 
Clearly, for a large number of modes, this procedure 
becomes impractical. 

Improved ROM Development Processes 

An outline of the improved simultaneous modal ex- 
citation ROM development process with the recent 
enhancements is as follows: 

1. Generate the number of functions (from a se- 
lected family) that corresponds to the number of struc- 
tural modeshapes; 

2. Apply the generated input functions simulta- 
neously via one CFD execution; these responses are 
computed directly from the restart of a steady rigid 
CFL3D solution (not about a particular dynamic pres- 
sure); 

3. Using the simultaneous input/output responses, 
identify the individual impulse responses using the 
PULSE algorithm (within SOCIT); 

4. Transform the individual impulse responses gen- 
erated in Step 3 into an unsteady aerodynamic state- 
space system using the ERA (within SOCIT); 

5. Evaluate/ validate the state-space models gener- 
ated in Step 4 via comparison with CFD results (i.e. , 
ROM results vs. full CFD solution results); 

An important difference between the original ROM 
process and the improved ROM process is stated in 
steps (2) of the outlines above. For the original ROM 
process, if a static aeroelastic condition existed, then 
a ROM was generated about a selected static aeroe- 
lastic condition. So a static aeroelastic condition of 
interest was defined (typically a dynamic pressure) and 
that static aeroelastic condition was computed using 
CFL3D as a restart from a converged steady, rigid so- 
lution. Once a converged static aeroelastic solution 
was obtained, the ROM process was applied about 
that condition. This implies that the resultant ROM 
is, of course, limited in some sense to the neighbor- 
hood of that static aeroelastic condition. Moving ’’too 
far away” from that condition could result in loss of 
accuracy. 

The reason for generating ROMs in this fashion 


was because no method had been defined to enable 
the computation of a static aeroelastic solution using 
a ROM. Any ROMs generated in this fashion were, 
therefore, limited to the prediction of dynamic re- 
sponses about a static aeroelastic solution including 
the methods by Raveh 11 and by Kim et al. 10 The im- 
proved ROM method, however, includes a method for 
generating a ROM directly from a steady, rigid solu- 
tion. As a result, these improved ROMs can then be 
used to predict both static aeroelastic and dynamic 
solutions for any dynamic pressure. In order to cap- 
ture a specific range of aeroelastic effects (previously 
obtained by selecting a particular dynamic pressure), 
the improved ROM method relies on the excitation 
amplitude to excite aeroelastic effects of interest. The 
details of the method for using a ROM for computing 
both static aeroelastic and dynamic solutions is pre- 
sented in another reference by the present author. 26 
For the present results, all responses were computed 
from the restart of a steady, rigid CFL3D solution, 
bypassing the need (and additional computational ex- 
pense) to execute a static aeroelastic solution using 
CFL3D. 

Simultaneous Excitation Input Functions 

In the situation where the goal is the simulta- 
neous excitation of a multiple-input multiple-output 
(MIMO) system, system identification techniques 34-36 
dictate that the nature of the input functions used 
to excite the system must be properly defined if ac- 
curate input /output models of the system are to be 
generated. The most important point to keep in mind 
when defining these input functions is that these func- 
tions need to be different, in some sense, from each 
other. This makes sense since, if the excitation inputs 
are identical and they are applied simultaneously, it 
becomes practically impossible for any system iden- 
tification algorithm to relate the effects of one input 
on a given output. This, in turn, makes it practically 
impossible for that algorithm to extract the individ- 
ual impulse responses for each input/output pair. As 
has already been well established, the individual im- 
pulse responses for each input /output pair are nec- 
essary ingredients towards the development of state- 
space models. With respect to unsteady aerodynamic 
MIMO systems, these individual impulse responses 
correspond to time-domain generalized aerodynamic 
forces (GAFs), critical to understanding unsteady 
aerodynamic behavior. The Fourier-transformed ver- 
sion of these GAFs are the frequency-domain GAFs 
which provide an important link to more traditional 
frequency-domain-based aeroelastic analyses. 

The question is how different should these input 
functions be and how can we quantify a level of ’’dif- 
ference” between each input function? Kim 10 uses the 
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familiar step function as the input function but with 
each step input being applied at a different point in 
time (lagged) in order to maintain some difference be- 
tween the input signals. Kim refers to this family of 
lagged step functions as ” almost orthogonal” . As Kim 
points out, the greater the lag between successive step 
inputs, the greater the difference between input sig- 
nals. This then presents a spectrum of possibilities 
where, at one end, one has step input functions with 
no lags (identical functions, not orthogonal) and at 
the other end the step functions are separated by an 
extremely large lag (in the limit, these are the indi- 
vidual step inputs from the original ROM process). 
Raveh 11 uses filtered white Gaussian noise signals as 
the input functions. In order to maintain the neces- 
sary level of difference amongst each input function, 
these functions need to be i.i.d. (independent, iden- 
tically distributed). The application of i.i.d. ran- 
dom functions for identification of MIMO systems is a 
traditional approach within the system identification 
community, especially in experimental settings. Since 
orthogonality (linear independence) is the most pre- 
cise mathematical method for guaranteeing the ’’dif- 
ference” between signals, the present research focuses 
on the application of families of orthogonal functions 
as candidate input functions. Using orthogonal func- 
tions directly provides a mathematical guarantee that 
the input functions are as different as mathematically 
possible. These orthogonal input functions then can 
be considered optimal input functions for the identifi- 
cation of a MIMO system. 

The four families of functions to be applied towards 
the efficient identification of a CFD-based unsteady 
aerodynamic state-space model are: lagged step, block 
pulse, Haar, and Walsh. Although it has been stated 
by Kim that the lagged-step functions are not com- 
pletely orthogonal (almost orthogonal), these func- 
tions are used in the present study for comparison 
purposes. The other three functions each represent 
an orthogonal family of functions. The block pulse 
functions, presented in Figure 3, are orthogonal and 
resemble a modified step input. An advantage of step- 
like functions is the broad frequency bandwidth that is 
excited by the impulsive nature of these functions. An 
added benefit of the block pulse functions over tradi- 
tional step functions is that block pulse inputs provide 
excitation to the system in both positive and nega- 
tive directions. From an unsteady aerodynamics point 
of view, it is important that modes of realistic wing 
configurations (e.g., with a non-symmetric airfoil) be 
excited in both directions for a more complete capture 
of relevant dynamics. 

The Haar functions, presented in Figure 4, com- 
prise one of the simplest form of wavelet functions. 
These functions are orthogonal by the mathematical 



t 


Fig. 3 Block pulse functions. 

definition of wavelet functions. Likewise, this family 
of functions has a similarity to step inputs and there- 
fore embodies the impulsive (i.e. , beneficial) nature of 
step inputs with regards to frequency bandwidth. The 
final set of orthogonal functions is the Walsh function, 
presented in Figure 5. 



t 

Fig. 4 Haar (wavelet) functions. 

It is important to mention that the Haar functions 
are generated using formulas based on a power of 
two algorithm. These functions have yielded excellent 
results when applied towards the identification of a 
CFD-based unsteady aerodynamic state-space model 
using a single CFL3D execution for a multi-mode con- 
figuration. However, due to the fact that these func- 
tions are generated based on a power of two algorithm, 
the number of modes drives the size and, subsequently, 
the record length of the input functions that are then 
input to CFL3D. As a result, this can severely limit the 
number of modes that can be excited at one time. For 
example, if ten modes are to be excited simultaneously, 
then the record length of the ten Haar functions gener- 
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Fig. 5 Walsh functions. 


ated (one per mode) is defined by the record length of 
the Haar function that corresponds to mode ten. Then 
taking two to the power of ten yields 1024 time steps. 
Twelve modes yields a record length for the Haar func- 
tions of 4096 time steps, thirteen modes yields 8192 
time steps, fourteen modes yields 16384 time steps, 
and so on. Clearly, for a very large number of modes, 
this family of functions is not practical. Therefore, 
only those results for the lagged step (not orthogonal), 
the block pulse (orthogonal), and the Walsh (orthog- 
onal) functions are presented. 

Supersonic Semi-Span CFD Model 

The configuration used for the present analyses is 
a supersonic semi-span wind-tunnel model known as 
the Rigid Semi-span Model (RSM) . This configuration 
has been tested several times at the NASA Lang- 
ley’s Transonic Dynamics Tunnel (TDT). The actual 
wind-tunnel model was fabricated using graphite and 
is very rigid. However, as part of collaborative stud- 
ies between the NASA Langley Research Center and 
the Boeing Company (Seattle), a ’’softened” computa- 
tional model of the RSM was developed by Dr. Moeljo 
Hong (Boeing). 7 The model was ’’softened” by simply 
reducing the four modal frequencies by a factor of four. 

The results presented are for Euler (invis- 
cid) solutions at a Mach number of 0.7 and an angle of 
attack of 3 degrees. This configuration does not have a 
symmetric airfoil and will, therefore, generate a static 
aeroelastic response in addition to the dynamic aeroe- 
lastic response. The method for using an unsteady 
aerodynamic state-space ROM for static aeroelastic re- 
sponses is not presented in this paper but is presented 
in another paper by this author. 26 Figure 6 presents 
the surface grid for the CFL3D RSM configuration. 

Results 

As previously stated, when using the lagged-step 
functions as excitation inputs, it is important to con- 
sider the lag (number of time steps) between each 
successive step input. It is the time lag between each 



Fig. 6 Computational grid of supersonic semi-span 
configuration. 

step input that determines the quality of the resultant 
model that is identified. A time lag of fifteen (15) time 
steps was used initially in order to explore the effect of 
the time lag on the resultant ROM. Presented in Fig- 
ure 7 is the CFL3D response in the first generalized 
coordinate due to the excitation of the four lagged- 
step inputs. The effect of each of the lagged steps is 
evident in the response of the generalized coordinate 
as evidenced by the four pulses in the response. 



Fig. 7 CFL3D response in the first generalized 
coordinate due to excitation using four lagged-step 
inputs with a lag of 15 time steps between each 
input. 
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Presented in Figure 8 is the CFL3D response in the 
first generalized coordinate due to the excitation of 
four lagged-step inputs but this time with a relative 
time lag of 50 time steps. The effect of each of the 
lagged steps is evident in the response of the general- 
ized coordinate as evidenced by the four pulses in the 
response. 



Fig. 8 CFL3D response in the first generalized 
coordinate due to excitation using four lagged-step 
inputs with a lag of 50 time steps between each 
input. 

Figure 9 presents the four block pulse functions 
input to the CFL3D RSM model for simultaneous ex- 
citation of the four modes. 

Presented in Figure 10 is the CFL3D response in 
the first generalized coordinate due to the excitation of 
the block pulse inputs. The effect of each of the block 
pulse inputs is evident in the response of the general- 
ized coordinate as evidenced by the four pulses in the 
response. Note that the block pulse functions were de- 
fined over an interval of 2000 time steps. The lagged 
step inputs were defined depending on the desired lag 
between each step input but all solutions (lagged step, 
block pulse, and Walsh) were executed for the same 
number of time steps. 

In Figure 11, Figure 12, Figure 13, and Figure 14 
each of the four Walsh functions is presented sepa- 
rately for clarity. The CFL3D response in the first 



Fig. 9 Block pulse functions input to the CFL3D 
RSM model for modal excitation. 

generalized coordinate due to the Walsh function in- 
puts is presented in Figure 15. Like the block pulse 
functions, the Walsh functions were generated over a 
record length of 2000 time steps. 

Processing the outputs for all four generalized coor- 
dinates for each input function through the PULSE al- 
gorithm in SOCIT, the time-domain generalized aero- 
dynamic forces (GAFs) can be identified. The first 
twenty-five time steps of the time-domain GAF in 
mode 1 due to an input in mode 1 (GAF 1,1) for the 
four input functions (lagged step (15 steps), lagged 
step (50 steps), block pulse, and Walsh) is presented 
as Figure 16. Clearly, if the focus is on the initial por- 
tion of the GAF (high frequency), then all four input 
functions yield identical GAFs. The results for the 
other GAFs (fifteen additional GAFs) is similar to the 
one shown in Figure 16. 

However, an analysis of the time-domain GAF 1,1 
from the 20th to the 200th time step, presented in Fig- 
ure 17, reveals noticeable differences due to the four 
input functions. Note that the vertical scale has been 
greatly expanded from that of Figure 16. This vari- 
ance can be considered an indication of the difficulty 
that the system identification algorithm is encounter- 
ing in attempting to identify the system. It can be 
seen that the greatest variance is for the lagged step 
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Fig. 12 Walsh input function for the first and 
second modes. 


Fig. 10 CFL3D response in the first generalized 
coordinate due to excitation using block pulse func- 
tions. 




Fig. 11 Walsh input function for the first mode 
(shown separately). 


Fig. 13 Walsh input function for the first, second, 
and third modes. 
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Fig. 14 Walsh input function for the first, second, 
third, and fourth modes. 



Fig. 16 Initial portion of the time-domain GAF 
in mode 1 due to an input in mode 1 (GAF 1,1). 

input with 15 time step delays, compared to the oth- 
ers. It is interesting to note the similarity between the 
lagged step (50 time steps) and the block pulse input 
functions. 




Fig. 17 Variation of GAF 1,1 for the four input 

Fig. 15 CFL3D response in the first generalized functions from the 20th to the 200th time step, 
coordinate due to excitation using Walsh functions. 

Additional comparison of these results can be ob- 
tained via the analysis of the frequency-domain version 
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of the GAFs. The frequency-domain GAFs from each 
of four input functions are compared to the frequency- 
domain GAFs from the use of serial (one-at-a-time) 
step inputs. Comparison of the simultaneous input re- 
sponses with those from the one-mode-at-a-time does 
not imply that the one-mode-at-a-time results are the 
correct results. That is, it is quite possible that (and 
most probable) that the responses obtained from the 
simultaneous excitation will be different from the re- 
sponses obtained from the one-mode-at-a-time excita- 
tions due to potential nonlinear modal coupling effects. 
The reason for comparing these two sets of responses is 
simply due to the fact that the one-mode-at-a-time re- 
sponses represent the original ROM approach. Appli- 
cation of an FFT to each one of the time-domain GAFs 
directly yields the desired frequency-domain GAFs. 8 
Only a subset of the sixteen GAFs will be presented in 
the figures below. The first comparison is for the real 
component of GAF 1,1, due to the lagged step input 
with 15 time-step delays, presented in Figure 18. The 
effect of the variance exhibited in the time domain is 
manifested as a difference at low reduced frequencies. 



Fig. 19 Imaginary component of GAF 1,1 for the 
lagged step (15 time steps) and serial step inputs. 



Fig. 18 Real component of GAF 1,1 for the lagged 
step (15 time steps) and serial step inputs. 


variations from the serial inputs are noticed as well, 
in particular with respect to the steady (zero reduced 
frequency) value. 



The imaginary component of GAF 1,1 due to the 
lagged step with 15 time-step delays is presented in 
Figure 19. The comparison for the imaginary com- 
ponent is quite good, with slight variations at certain 
frequency ranges. 

The real and imaginary components of GAF 1,1 due 
to the lagged step with 50 time-step delays are pre- 
sented in Figure 20 and Figure 21. Here again, some 


Fig. 20 Real component of GAF 1,1 for the lagged 
step (50 time steps) and serial step inputs. 

The real and imaginary components of GAF 1,1 due 
to the block pulse function inputs are presented as Fig- 
ure 22 and Figure 23. For both the real and imaginary 
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Fig. 21 Imaginary component of GAF 1,1 for the 
lagged step (50 time steps) and serial step inputs. 

components, noticeable improvement in the compari- 
son is achieved with the use of the block pulse function. 


Fig. 23 Imaginary component of GAF 1,1 for the 
block pulse and serial step inputs. 

comparison is achieved with the use of the Walsh func- 
tion. In particular, for both the block pulse and the 
Walsh functions, the steady value (value at zero re- 
duced frequency) is well captured. 



Fig. 22 Real component of GAF 1,1 for the block 
pulse and serial step inputs. 

The real and imaginary components of GAF 1,1 due 
to the Walsh function inputs are presented as Figure 24 
and Figure 25. For both the real and imaginary com- 
ponents, once again, noticeable improvement in the 



Fig. 24 Real component of GAF 1,1 for the Walsh 
and serial step inputs. 

A comparison of the frequency-domain GAF 2,2 is 
presented for the lagged step (50 time-step delays) and 
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Fig. 25 Imaginary component of GAF 1,1 for the 
Walsh and serial step inputs. 

the Walsh function inputs. Figure 26 and Figure 27 
contain the comparison of the real and imaginary com- 
ponents of GAF 2,2 for the lagged step (50 time-step 
delays) and the serial (step) inputs. Some differences 
are noticed, in particular in the low reduced frequency 
range for both the real and imaginary components. 



Fig. 26 Real component of GAF 2,2 for the lagged 
step (50 time-step delays) and serial step inputs. 


Fig. 27 Imaginary component of GAF 2,2 for the 
lagged step (50 time-step delays) and serial step 
inputs. 

Figure 28 and Figure 29 are the comparison of the 
real and imaginary components of GAF 2,2 for the 
Walsh and the serial (step) inputs. An improvement 
in the comparison is once again noted. 

The time-domain GAFs are then processed through 
the ERA (part of SOCIT) in order to generate a state- 
space model of the unsteady aerodynamic system. The 
resultant state-space models (for the lagged step with 
50 time-step delays and the Walsh functions) were 12th 
order in dimension. To demonstrate the effectiveness 
of the new input functions, a representative result is 
presented as Figure 30, a comparison of responses in 
the four generalized coordinates for the CFL3D solu- 
tion and the ROM solution. The CFL3D responses are 
from a CFL3D solution that includes both static and 
dynamic aeroelastic responses simultaneously. Typ- 
ically, static and dynamic aeroelastic responses are 
computed separately in order to reduce computational 
cost. A static aeroelatic solution is computed as the 
restart solution from a steady, rigid solution using an 
artifically high structural damping value in order to 
accelerate convergence to a static aeroelastic solution. 
Then, a dynamic aeroelastic solution is computed as 
the restart solution from the converged static aeroe- 
lastic solution. In the paper by the present author, 26 
a method is introduced that enables the prediction 
of static and dynamic aeroelastic solutions from the 
same ROM. For the present paper, the sample results 
presented in Figure 30 are for a full solution; that is, 
a solution that includes both static and dynamic re- 
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Fig. 28 Real component of GAF 2,2 for the Walsh 
and serial step inputs. 


sponses for a given dynamic pressure. In this case, the 
results presented are for a dynamic pressure of 0.1 psi. 
Notice that the solutions start at a generalized coordi- 
nate value of zero (static, rigid solution) and converge 
to a non-zero mean value (static aeroelastic solution) 
with the dynamic solution superimposed on the static 
solution. As can be seen, the comparison is very good 
and the ROM captures the combined static and dy- 
namic aeroelastic responses predicted by the CFL3D 
solution. 



, xKT 




1000 

Time Steps 



Fig. 29 Imaginary component of GAF 2,2 for the 
Walsh and serial step inputs. 


Fig. 30 Full solution (static plus dynamic aeroe- 
lastic responses) generalized coordinate responses 
at a dynamic pressure of 0.1 psi. 


Concluding Remarks 

An improved method for generating unsteady aero- 
dynamic state-space reduced-order models (ROMs) 
was introduced that enables the simultaneous excita- 
tion of an unsteady aerodynamic CFD model, con- 
sisting of any number of modes, with a single CFD 
execution. This is achieved by the introduction of or- 
thogonal functions as the input excitations to the sys- 
tem. Comparisons of generalized aerodynamic forces 
(GAFs) were presented for the lagged step input, block 
pulse input, and the Walsh function input. Results 
compared well with the results from the one-mode- 
at-a-time method that, for four modes, required four 
separate CFD executions. Additional research is un- 
derway to continue the evaluation of these functions 
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on more complex configurations with a larger number 
of modes in order to define the optimal input (and 
input amplitudes) for the simultaneous excitation of 
unsteady aerodynamic systems. In addition, research 
will be performed in the system identification arena in 
order to minimize the error associated with the iden- 
tified state-space models. 
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